Why the a\ meson might be a difficult messenger for the restoration of chiral 

symmetry 
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We perform a theoretical analysis of the a\ resonance mass spectrum in ultra-relativistic heavy 
ion collisions within a hadron/string transport approach. Predictions for the ai yield and its mass 
distribution are given for the GSI-FAIR and the critRHIC energy regime. The potential of the a\ 
meson as a signal for chiral symmetry restoration is explored. In view of the latest discussion, we 
investigate the decay channel a\ — > 771 in detail and find a strong bias towards low ai masses. This 
apparent mass shift of the ai if observed in the 77r channel might render a possible mass shift due 
to chiral symmetry restoration very difficult to disentangle from the decay kinematics. 
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One of the main goals of relativistic heavy ion physics 
is to reach densities and temperatures high enough to re- 
store chiral symmetry [l|. Chiral symmetry is a symme- 
try of quantum chromodynamics, which is exact if quark 
masses are zero and approximate if quark masses are 
small. It is spontaneously broken in nature, but expected 
to be restored at sufficiently high densities and temper- 
atures. The restoration of chiral symmetry implies a 
change in the spectral functions of vector mesons (e.g. 
the p meson) and leads to a degeneracy of the spectral 
functions of the p and its chiral partner, the a\ meson. 
This means that the masses of the chiral partners become 
equal in the case of full chiral symmetry restoration or 
approach each other in the case of a partial restoration 
of the symmetry. 

Especially the recently observed broadening of the p 
meson spectral function by the NA60 collaboration and 
the corresponding di-lepton mass spectrum has been 
interpreted as a signal of chiral symmetry restoration 
[lilS IH IE HI- In fact, the NA60 collaboration measured 
the p meson spectral function in In+In systems at the 
highest SPS energy of 158 AGeV and observed a devia- 
tion from the vacuum Breit-Wigner distribution . This 
has triggered various theoretical investigations [all, H(| • 
In summary these studies suggest that some in-medium 
effects have to be considered, but a conclusive interpre- 
tation of the data is still under discussion. 

Also the HADES collaboration has recently presented 
first results on di-electron spectra in light systems at low 
beam energies (C+C at 2 AGeV) [ll|. Here a deviation 
from the hadronic vacuum cocktail is visible in the mass 
region of 500 to 700 MeV. This has been discussed as a 
possible observation of partial chiral symmetry restora- 
tion and the resulting change in the p meson spectral 
function. Similar data has also been measured by the 
CERES experiment at CERN in massive nuclear reac- 
tions at high energy [l2| . In spite of the ongoing experi- 
mental and theoretical efforts, there are numerous effects 
that have to be taken into account for a full understand- 
ing of data. Thus, it is questionable that a mass shift of 
the p meson alone can be regarded as a "smoking gun" 



signal of chiral symmetry restoration [l3|, [l4j]. There- 
fore a more robust signature of chiral symmetry restora- 
tion is needed. Theory predicts that in the case of a 
full restoration of the symmetry the spectral functions of 
the p meson and its chiral partner the a\ meson become 
degenerate. The important point is that this statement 
is independent of any mass shift or broadening. Thus, 
it has been proposed to measure the a\ mass spectrum 
in a hot and dense medium and compare it to the mass 
spectrum of the p meson (l5| . If the degeneracy would 
be observed it is expected to serve as an unambigious 
experimental signal for the detection of chiral symmetry 
restoration in the hot and dense medium. 

In this paper we argue that the measurement of the a\ 
(1260) meson may not result in straightforward insights 
for the understanding and the detection of the chirally 
restored phase. We discuss the decay kinematics of the 
a\ meson and argue that an apparent mass shift, or re- 
spectively a broadening of the mass spectrum originates 
from mass dependent branching ratios. This effect is not 
unique to heavy ion reactions, but is nearly independent 
of energy and system size. Furthermore we predict a\ 
mass spectra for Au+Au and p+p collisions at 20 and 
30 AGeV beam energy. The p+p calculations can serve 
as a vacuum reference. These systems and energies are 
experimentally accessible at FAIR-GSI, NA61 and the 
critRHIC program in the near future. 

For our calculations we utilize the UrQMD(v2.3) 
model, a non-equilibrium transport approach, which re- 
lies on the covariant Boltzmann equation. All cross sec- 
tions are calculated by the principle of detailed balance 
and the additive quark model or are fitted to available 
data. UrQMD does not include any explicit in-medium 
modifications or effects to describe the restoration of chi- 
ral symmetry. The model allows to study the full space 
time evolution of all hadrons, resonances and their de- 
cay products in hadron-hadron or nucleus-nucleus colli- 
sions. This permits to explore the emission patterns of 
resonances in detail and to gain insight into their origins 
and decay channels. For previous studies of resonances 
within this model see [II HE H3, H GJ]. For further 
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details about the model the reader is referred to [2(| [2l| . 

Experimentally, the reconstruction of resonances is 
challenging. One usual technique is to reconstruct the 
invariant mass spectrum for single events. Then, an in- 
variant mass distribution of mixed events is generated 
(here, the particle pairs are uncorrelated by definition). 
This distribution is substracted from the invariant mass 
spectrum of the single (correlated) events. As a result 
one obtains the mass distributions and yields (after all 
experimental corrections) of the resonances by fitting the 
resulting distribution with a suitable function (usually a 
Breit-Wigner function peaked around the pole mass of 
the respective resonance) [H, [23[ . If a daughter parti- 
cle (re-) scatters before reaching the detector the signal 
for the experimental reconstruction is lost. Especially 
for strongly interacting decay products this effect can 
be sizeable. In addition, due to the statistical nature of 
the reconstruction, detailed information on the particle 
properties and their origin are difficult to obtain. Also 
possible deviations from a Breit-Wigner distribution can 
be overseen due to the dependence on the background 
substraction. 

Thus, we apply a different technique for the extraction of 
resonances from the model. We follow the indiviual decay 
products of each decaying resonance (the daughter par- 
ticles). If the daughter particles do not rescatter in the 
further evolution of the system, the resonance is counted 
as "reconstructable" . The advantage of this method is 
that it allows to trace back the origin of each individual 
resonance to study their spatial and temporal emission 
pattern. It also allows to explore the reconstruction effi- 
ciency in different decay branches. 

The decay channels of the a\ meson are not fully ex- 
perimentally investigated and details of the branching ra- 
tios are unknown [24| . However, a most promising decay 
channel for the investigation of the restoration of chiral 
symmetry seems to be the decay a\ — > 771", due to the fact 
that the photon does basically not interact with the sur- 
rounding (hadronic) medium. A study of all other decay 
channels would imply to study three particle correlations 
or respectively correlations between resonances and sta- 
ble particles, which is very tedious, if not impossible in 
large systems — however, see also the discussion at the 
end of this letter. Thus, experimentally, the a\ — > 77T 
channel seems the only feasible candidate to measure the 
Oi meson in heavy ion collisions. 

One problem is that the branching ratio into this cer- 
tain decay channel is still unknown. Further experimen- 
tal studies in elementary systems are in order to obtain 
more precise quantitative results in theoretical investiga- 
tions. While the observed branching ratio r ai ^ 77r /T^°* 
~ 10~ 3 [24[ we use an increased branching ratio of 10 _1 
in this letter to improve statistics. The shape of the spec- 
tra and the discussion of the decay kinematics are essen- 
tially independent of the exact value of the branching 
ratio however one should keep this upscaling in mind for 
the absolute multiplicities observable in the "fir channel. 

In Fig. [T] the mass spectrum of the a\ meson for 



p+p collisions at 20 and 30 AGeV obtained from the 
UrQMD calculation is shown. One observes a clear peak 
around the pole mass of the a\ meson. At lower masses 
a shoulder-like structure is visible, which will be of im- 
portance in the discussion of the di — > 77r channel later 
on. Note that this mass spectrum is narrower than the 
one obtained experimentally from r decays [25j | . see also 
26] for a detailed discussion of the a\ spectral shape. 
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FIG. 1: Mass distribution of ai mesons in proton-proton col- 
lisions at Ei a b = 20 and 30 AGeV as obtained from UrQMD 
calculations. 

In Fig. [2] the mass spectrum of a\ mesons for central 
(b < 3.4 fm) Au+Au collisions at 20 and 30 AGeV as 
obtained from UrQMD calculations is depicted. As in 
the p+p case one observes a peak around the pole mass 
(although slightly shifted to lower masses for kinematic 
reasons discussed in [l8j]) and a shoulder- like structure 
at lower masses of roughly 500 MeV. This shoulder orig- 
inates from decays into the jtt channel and will be dis- 
cussed in the following. 

What happens now if one explicitly triggers on the 
decay channel a\ — > 77T that seems most suitable for the 
study of the a\ in heavy ion reactions? As discussed in 
[2CI |27| the branching ratios of resonances depends on 
the mass of the decay products. The total decay width 
Ttot(M) of a resonance is defined as the sum of all partial 
decay widths and depends on the mass of the excited 
resonance: 

r tot (M) = r iAM), (i) 

br={i,j} 

where Tij(M) is the partial decay width, M is the mass 
of the resonance and the summation over Nb r denotes a 
summation over all possible decay channels. The partial 
decay widths Tij(M) for the decay into the exit channel 
with particles i and j is given by [2(3, H3| : 
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FIG. 2: Mass distribution of a± mesons in central Au+Au 
collisions (b < 3.4 fm) at 20 and 30 AGeV as obtained from 
UrQMD calculations. 
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here Mr denotes the pole mass of the resonance, T R J 
partial decay width into the channel i and j at the pole 
and I the decay angular momentum of the exit channel. 
{pij(M)} denotes the momentum of the decay products 
in the center of momentum frame. 

Fig. [3] shows the branching ratio of the a\ meson 
as a function of the mass of the a\ as obtained from 
UrQMD calculations, wherein the definitions Q] and [2] are 
implemented including the finite width of decay particles. 
Filled squares depict the branching ratio of a\ mesons 
into the exit channel /97r, whereas open squares depict 
the branching ratio into 777. Also shown is a normalized 
Breit-Wigner distribution (full line) and the normalized 
mass spectrum of the a\ meson as obtained from proton- 
proton collisions at 20 AGeV from UrQMD. One observes 
that at masses lower than 750 MeV the decay channel 
of a\ — y 77T dominates because the decay channel into 
pir is kinematically suppressed. At masses greater than 
750 MeV the pw decay channel is dominantly populated 
and the contribution from the 777 channel becomes less 
important. Depicted in Fig. [3] are only two of the possi- 
ble decay channels listed in [2J|. All other exit channels 
consist of even heavier decay products. 

Fig. 2] depicts the mass spectrum (with arbitrary nor- 
malization) as obtained from folding the Breit-Wigner 
distribution of the a\ meson with the mass dependent 
branching ratio as depicted in Fig. [3J One observes a 
double peaked structure, with one peak at the nominal 
pole mass of 1260 MeV and an additional peak at much 
lower masses. Also note that it is not completely indepen- 
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FIG. 3: Mass dependent branching ratios for the a± me- 
son with the two exit channel of ^tt and piv as calculated 
from UrQMD. Filled squares depict the branching ratio of en 
mesons into the exit channel pn, whereas open squares depict 
the branching ratio into 'yn. Below a mass of 750 MeV the 
decay channel a 1 — > pn is kinematically suppressed and the 
channel a± — > 771- dominates. At masses above 750 MeV the 
branching ratio into pn increases steeply. The grey shaded 
area depicts a normalized Breit-Wigner distribution around 
the ai pole mass, whereas the circles depict the normalized 
mass spectrum of the ai meson as obtained from UrQMD 
calculations for p+p collisions at 20 AGeV. 



dent of the partial decay widths, i.e. the branching ra- 
tios. In this particular case a total width of 400 MeV was 
assumed with the (enhanced) partial widths of 40 MeV 
(ai — -> 771") and 360 MeV (ai — > pw). 




FIG. 4: Sketch of the mass spectrum for the ai meson as 
obtained from folding the Breit-Wigner distribution with the 
branching ratios as shown in Fig. 3. One observes a double 
peak structure with one peak centered around the nominal 
pole mass of 1260 MeV and another peak at lower masses. 

After these semi-quantitative discussions, it is clear 
that a non-trivial a\ mass spectrum has to be expected 
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in the full UrQMD calculation, if a trigger on the jit exit 
channel is employed. Let us therefore check this effect 
within the full transport model calculation. Fig. [5] de- 
picts the mass spectrum for those a\ mesons which can 
be reconstructed in the aj — > jtt decay channel. One 
observes a clear double peak structure, with one peak 
centered around the pole mass and one peak in the range 
of 400-600 MeV. The low mass peak is even more promi- 
nent in p+p than in Au+Au collisions. Thus, a possible 
ai mass shift due to chiral symmetry restoration might 
be difficult to distinguish from a scenario without mass 
shift but including mass dependent branching ratios. 
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FIG. 5: Mass distribution of a\ mesons which can be recon- 
structed in yir correlations in nucleus-nucleus and proton- 
proton collisions at 20 and 30 AGeV. Note that the p+p 
curves have been multiplied by 100 for visibility, and all cal- 
culations use the enhanced branching ratio into 771". 



Measuring the a\ from a correlations of p mesons 
and pions might be more robust, however is experimen- 
tally even more demanding. The pir channel will avoid 
the problem of mass dependent branching ratios. In- 
deed, the HADES (and later the CBM) experiments offer 
the unique possibility to measure correlations between p 
mesons and pions, where the p mesons are reconstructed 
via the decay channel p° — > e + e— . 

These measurements might indeed provide a novel and 
up to now unexplored route to obtain insights into the 
transition from the chirally broken to the chirally re- 
stored phase. 

In summary, we have presented a\ mass distributions 
for p+p and Au+Au reactions at Ei a b = 20 and 30 AGeV. 
It was shown that the reconstruction of the a\ mass dis- 
tribution in the jit channel is strongly biased towards low 
masses. This can be traced back to the strong mass de- 
pence of the a\ branching ratios. Thus, we conclude that 
a possible observation — e.g. at the planned GSI-FAIR 
facility or in a future critRHIC run — of a modified a\ 
mass distribution measured in the a\ —* jit exit channel 
cannot unambigiously signal an approach towards chiral 
symmetry restoration. However a measurement in the 
di — » pit channel as it is possible with the HADES and 
CBM experiments might lead to new information about 
the approach to chiral symmetry restoration. 

The computational resources have been provided by 
the Center for Scientific Computing, CSC at Frankfurt 
University. We thank Volker Koch for fruitful discus- 
sions. This work has been financially supported by GSI 
and BMBF. S.V. thanks the Helmholtz foundation for 
additional financial support. 



[1] S. A. Bass, M. Gyulassy, H. Stoecker and W. Greiner, J. [14] 
Phys. G 25, Rl (1999) 

V. Koch and G. E. Brown, Nucl. Phys. A 560, 345 (1993). [15 

G. E. Brown and M. Rho, Phys. Rept. 269, 333 (1996) [16 
W. Cassing, E. L. Bratkovskaya, R. Rapp and [17 
J. Wambach, Phys. Rev. C 57, 916 (1998) 

R. Rapp and J. Wambach, Adv. Nucl. Phys. 25 1 (2000) [18 
W. Cassing and E. L. Bratkovskaya, Phys. Rept. 308, 65 

(1999) [19 

R. Arnaldi et al. [NA60 Collaboration], Phys. Rev. Lett. [20 

96, 162302 (2006) [21 

J. Ruppert, T. Renk and B. Muller, Phys. Rev. C 73, [22 
034907 (2006) 

J. Ruppert and T. Renk, Eur. Phys. J. C 49, 219 (2007) [23 

H. van Hees and R. Rapp, Phys. Rev. Lett. 97, 102301 [24 
(2006) 

G. Agakishiev [HADES Collaboration], Phys. Rev. Lett. [25 
98, 052302 (2007) 

D. Miskowiec [CERES Collaboration], Nucl. Phys. A [26 

774, 43 (2006) [27 
S. Vogel and M. Bleicher, Phys. Rev. C 74, 014902 (2006) 



D. Schumacher, S. Vogel and M. Bleicher, Acta Phys. 
Hung. A 27, 451 (2006) 

Volker Koch, workshop on dileptons at CBM, GSI, 2007 
M. Bleicher and J. Aichelin, Phys. Lett. B 530, 81 (2002) 
M. Bleicher, Nucl. Phys. A 715, 85 (2003) 
[arXiv:hep-ph/0212378] . 

M. Bleicher and H. Stoe cker, J. Phys. G 30, Sill (2004) 
arXiv:hep-ph/0312278j . 



Vogel and M. Bleicher, arXi v:nucl-th/ 0505027 
S. A. Bass et al, Prog. Part. Nucl. Phys. 41, 225 (1998) 
M. Bleicher et al, J. Phys. G 25 1859 (1999) 
J. Adams et al. [STAR Collaboration], Phys. Rev. C 71, 

064902 ( 2005) 

R. Witt, |arXiv:nucl-ex/0701063| 

W. M. Yao et al. [Particle Data Group], J. Phys. G 33, 
1 (2006). 

S. Schael et al. [ALEPH Collaboration], Phys. Rept. 421, 

191 (2005) 

M. Wagner and S. Leupold. larXiv:0708.2223l [hep-ph] . 
H. Sorge, Phys. Rev. C 52, 3291 (1995) 



